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Abstract 

A model for the twist-3 wave function 'iljp{x,]<.±) of the pion has been constructed based on 
the moment calculation by applying the QCD sum rules, whose distribution amplitude has a 
better end-point behavior than that of the asymptotic one. With this model wave function, the 
twist-3 contributions including both the usual helicity components (Ai -|- A2 = 0) and the higher 
helicity components (Ai -|- A2 = ±1) to the pion form factor have been studied within the modified 
pQCD approach. Our results show that the twist-3 contribution drops fast and it becomes less 
than the twist-2 contribution at ~ lOGeV^. The higher helicity components in the twist-3 
wave function will give an extra suppression to the pion form factor. The model dependence of 
the twist-3 contribution to the pion form factor has been studied by comparing four different 
models. When all the power contributions, which include higher order in q^, higher helicities, 
higher twists in DA and etc., have been taken into account, it is expected that the hard con- 
tributions will fit the present experimental data well at the energy region where pQCD is applicable. 
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I. INTRODUCTION 



The most challenging problems for applying the perturbative QCD (pQCD) to exclusive 
processes have long been discussed and analyzed in many papers, such as the pQCD ap- 
plicability to the exclusive processes at experimentally accessible energy region due to the 
end-point singularity; to estimate the contributions from power corrections, which includes 
higher order in a^, higher helicities, higher twists in distribution amplitude (DA), higher 
Fock states and etc.; to estimate the uncertainties from perturbatively incalculable DAs. 

The pion form factor can be obtained through the definition 

{Ap)\JMp)) = ip + p),FAQ'), (1) 

where = J^i^-iqi'j^qi, with the quark flavor i and the relevant electric charge Cj, is the 
vector current. The momentum transfer = —Q^ = {p — p'Y is restricted in the space-like 
region. The pQCD applicability to the pion form factor at the experimentally accessible 
energy region has been raised by Ref. Q| and attracted much attention for many years. In 
the modified pQCD approach that is proposed in Ref. [^, i.e. the transverse momentum 
dependence (fc-r dependence) as well as the Sudakov corrections are taken into account in 
the calculations, we have the following factorization formula P, 0, 0, 0], 

n,m 

(2) 

where [dxid'kj\n is the relativistic measure within the n-particle sector, n, m extend over the 
low momentum states only and Tnm are the partonic matrix elements of the effective current 
operator. Here the helicity states of the pion are implied in both sides. The dependence on 
the scale separating low (non-perturbative) and high momenta (perturbative) is indicated 
by /ij. For the valence quark state of the pion, its light cone (LC) wave functions are defined 
in terms of the bilocal operator matrix element |6|], 



(vr(p)|g/3(^)ga|0) = ^ f^dx J d'k^e'^'P-^-''--^-^ 



X ||^75^.(x, k^) - ;x.75 (^p(x, k^) - a^.p^;^'- '^^^^^ j | , (3) 

where /i^ = m^/irriu + rrij) and is the pion decay constant, whose experimental value is 
130.7 ± 0.1 ± 0.36Mel^y . ij^{x, k_L) is the leading twist (twist-2) wave function, ^p{x, kj_) 



and ipa{x, k±) are sub-leading twist (twist-3) wave functions that correspond to the pseudo 



scalar structure and the pseudo-tensor structure respectively 
and the wave function ip{x, k±) are related by 



(f){x) 



-^(x,k_L). 



The distribution amplitude 



(4) 



It has been shown in different approaches]^, 0] that applying pQCD to the pion form 
factor begins to be self-consistent for a momentum transfer at about ~ 4GeV^. The 
next-to-leading order (NLO) QCD corrections to the pion form factor at large momentum 



transfer has also been analyzed 



Q,liilQ, 



13 



15 



16 



17l |. Ref.jl7| presents a complete 



NLO pQCD prediction for the pion form factor and it shows that a reliable pQCD prediction 
can be made at a momentum transfer around (5 — 10)GeV with corrections to the LO 
results being up to ~ 30%. The theoretical uncertainty related to the renormalization scale 
ambiguity has been estimated to be less than 10% and for all the considered DAs, concerning 
the choices of the renormalization schemes and the factorization scales, the ratio of the NLO 
to the LO contribution to the pion form factor F-j^lQ"^) is greater than 30% as < 20GeV^. 

A detailed calculation about the higher helicity components' contributions to the hard 
part and the soft part of the pion form factor within the LC pQCD approach was presented 
in Ref . . Their results show that by fully keeping the transverse momentum dependence 
in the hard part, the asymptotic behavior of the hard scattering amplitude from the higher 
helicity components is of order but it can give a sizable contribution to the pion form 

factor at the present experimentally accessible energy region. 

Other power corrections are from the higher twist structures in the pion DA. In the 
literature, based on the asymptotic behavior of the twist-3 DAs, especially (t>p^{x) = 1, most 



of calculations give large twist-3 contributions 



m 



3, 



21 



22 



23[, i.e. the twist-3 contribution 



to the pion form factor is comparable or even larger than that of the leading twist in a wide 
intermediate energy region, e.g. ~ (2 — AO)GeV^. It is hard to believe these results 
are reliable, since the power suppressed corrections make such a large contribution up to 
40GeV^. However, because the end-point singularity becomes more serious, the calculations 
for these higher twist contributions have more uncertainty than that for the leading twist. 
In fact, one may find that such kind of large contribution comes mainly from the end-point 
region and is model dependent. It means that one should try to look for a reasonable twist-3 
wave function with a better behavior in the end-point region than that of the asymptotic 
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one, and the twist-3 contribution might be less important and less uncertainty. 



Recently in Ref. 



24| , based on the moment calculation, the authors obtained a new form 



for which has a better behavior at the end-point re gion than that of the asymptotic 



one. Their approach is different from that of Refs.j8|, |25|, |26|, i.e. they did not apply the 
equation of motion for the quarks in the hadron and determined the coefficients of the 
Gegenbauer polynomial expansions directly from the DA moments obtained in the QCD 
sum rules. The 0p(x) obtained in Ref.j^l can be used to suppress the end-point singularity 
coming from the hard scattering kernel. In this paper, we will develop it to construct a 
model wave function ippi^x^k^) and apply it to calculate the twist-3 contributions to the 
pion form factor. 

The remainder of the paper is organized as follows. In Sec. II, we construct a model 
for the pionic twist-3 wave function k^) with the help of the moment calculation in 



Ref. 



2J|. And in Sec. Ill, the twist-3 contribution to the pion form factor, including those 



coming from the higher helicity components, will be studied within the modified pQCD 
approach. In Sec. IV, we discuss the model dependence for the twist-3 contribution. Finally 
we summarize our results and give the combined hard contributions to the pion form factor 
in Sec.V. 



II. A MODEL FOR THE PIONIC TWIST-3 WAVE FUNCTION 

For the twist-3 DAs, since the asymptotic behavior of 4>p{x) and (pa{x) are, 0p*(a;) ~ 1 
and ~ Qx{l — x) respectively, one may observe that the end-point singularity comes 

more seriously from (j)p{x) than from (f)cr{x)- With (f)cr{x) in the asymptotic form, the end- 
point singularity coming from the hard scattering kernel can be cured, while the asymptotic 
behavior of (f)p{x) can not suppress such kind of end- poin t singularity. 

The pion twist-3 DAs have been studied in Refs.0,|2^,0|- They employed the conformal 
symmetry and the equations of motion of the on-shell quarks within the hadron to get the 
relations among the two-particle twist-3 DAs, i.e. 0p(O 0cr(O (here and hereafter 
= {2x — 1)), and the three-particle twist-3 DA 03,r(tti) (o^i {i = 1,2,3) is the longitudinal 
momentum fraction of the corresponding constituent in the three-particle state (higher Fock 
state, e.g. \udg)) of the pion and satisfies Y^i^^i = !)• Then they took the moments of 
037r(ai) to obtain the approximate forms for the two-particle twist-3 DAs. However as has 
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been argued in Ref.j^l, since the quarks are not on-shell, it is questionable to use the 
equation of motion. So Ref. Q] suggested to calculate the moments of the pion two-particle 
twist-3 DAs directly from the QCD sum rules. 

Under the approximation that the lowest pole dominate and the higher dimension con- 
densates are negligible, the sum rule for the moments of 0p(O can be written as [24]. 



(5) 



where M is the Borel parameter and (^p") is the moment of (ppiC,), which is defined by 
(^p") = 1 /-I i'^^4>p{i)di- The parameter St, in Eq.® should be chosen to make the moments 
and the parameter most stable against M ^ in a certain range. In Eq. Q , one may observe 
that the usual /i7r-dependence in the sum rule for the moments of 0p(O|Sl ISSl |26| has been 
replaced by an undetermined parameter uiq. With the help of Eq.ijS)), setting {S,^) = 1 and 
varying the Borel parameter M in a reasonable range, we can obtain the values for the 
moments that are necessary to fit the parameters for our model wave function. 

Now we construct a model wave function ippi^^ of the twist-3 part that is related to 
(/)p(x) by the definition Eq.(j3]). The intrinsic transverse momentum dependence is determined 
by the non-perturbative dynamics and at present we cannot solve it. Ref.j2^ suggested a 
connection between the equal-time wave function ipc.m.{^-i-) ^^e rest frame and the LC 
wave function ipLci^, k±) in the infinite momentum frame, i.e. 

which expressed that the LC wave function should be a function of the bound state off-shell 



energy. Eq.Q is the so called BHL prescription [28|. Recently, some improvements on the 
transverse momentum dependence of the wave function have been given in Ref. [3], which 
presents a systematic study of the B meson LC wave function in the heavy-quark limit and 
by applying the QCD equations of motion. Their results show that under the Wandzura- 
Wilczek approximation 0, [3^ , the transverse and the longitudinal momenta in the B meson 
wave function are correlated through the combination ~ k5^/x(l — x)^. By adopting the 



^ Here x = uj/ (2A) g (0, 1), where lu, roughly speaking, is the longitudinal momentum of the light quark in 
B meson and A = {M — nib) is the "effective mass" of B meson in the heavy quark effective theory. 
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above prescription Eq. © and by using the harmonic oscillator model in the rest frame, the 
transverse momentum dependence part, S(a;, k^), can be written asjs^, 

where m and (3 are the quark mass and the harmonic parameter, respectively. Combining it 



with the new form of (f)p{C,), which is in the Gegenbauer polynomial expansion 
one can construct a model wave function with fc^ dependence. 



a, 



25 



M.M = (1 + BA'\^ - 2.) H- C,Cl'\, - 2.-))j^exp (-^^^i^) , (8) 

where 6*2^^(1 — 2x) and Cy^{l — 2x) are Gegenbauer polynomials and the coefficients Ap, 
Bp and Cp can be determined by the DA moments. In Eq.(jHI), only the first three terms 
in the Gegenbauer polynomial expansions have been considered. Since the higher moments 
of 0p(O obtained from the sum rule (Eq.(jS))) depends heavily on the Borel parameters, it 
is unreliable to do further expansions, so we only take the first three moments which have 
a better confidence level for our discussion. The parameters m and f3 can be taken from 
assuming the same kx dependence as the twist-2 wave function, and here we takej^^ 

m = 29QMeV, p = 385MeV, (9) 

which are derived for (k_|_^) ^ (356MeV^)^. From the model wave function Eq.(jHJ, we obtain 

MO = + BpC',/\0 + CpCy'iO) exp (- ^^2(7la ) ■ 

Reasonable ranges for the 0p(O nioments have been given in Ref. |i24] by applying the QCD 
sum rules (Eq.®), i.e. (^^^ _ (0.340,0.360) and (^^) ~ (0.160,0.210). Here we take 

(0 = 1, (e') = 0.350, (e') = 0.185, (11) 

for our latter discussion. The parameters in the wave function can then be determined as, 

Ap = 2.841 X 10-^Me1/-^ Bp = 1.302, Cp = 0.126. (12) 

As is shown in Fig.C), the shape of the present DA for 4'p{^) is very close to the one that 
is proposed in Ref.j24|. 

In the model wave function defined in Eq.(jH)), only the usual helicity components (Ai + 
A2 = 0) have been taken into account, while the higher helicity components (Ai + A2 = ±1) 
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FIG. 1: Different type of twist-3 DA. 



the asymptotic DA, the DAs of Ref. 
dash-dot line respectively. 



he solid line is for our 4>p{(,). And for comparison, we list 
24l | and Refs.P, 13] in diamond line, the dashed line and the 



which come from the spin-space Wigner rotation have not been considered. As has been 



pointed out in Refs. 



32| |. there is a large suppression coming from the higher hehcity 



components in the leading twist wave function, and one may expect that the higher helicity 
components in the higher twist wave functions also will do some contributions to the pion 
form factor. So we need to consider the higher helicity components in the twist-3 wave 
function. The full form for the LC wave function, i.e. ip^{x,\<.j_), which includes all the 
helicity components, can be found in the appendix. From ^^(x, k^), one may directly 
find that its DA 0^(0 is almost coincide with 0p(O and for simplicity, we can take the 
approximate relation, 0^(0 ~ 

In Fig.(^, we show our - 
totic DA, the DAs of Rel^ and Refs. 
dash-dot line, respectively. One may observe that the possible end-point singularity coming 
from the hard scattering kernel will be suppressed in our DA and the twist-3 contribution 
can be greatly suppressed at the present experimentally accessible energy region. 



[C,) in solid line, and for comparison, we also present the asymp- 
• ^, in the diamond line, the dashed line and the 
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III. THE TWIST-3 CONTRIBUTION TO THE PION FORM FACTOR IN THE 
MODIFIED PQCD APPROACH 



In the large region, by considering only the lowest valence quark state of the pion (i.e. 
n = m = 2 m Eq.(j2I)) and by doing the Fourier transformation of the wave function with 
the formula. 



?/'(xi,k_L;/i/) 



(2^ 



jb-k I 



V^(xi,b;/i/) 



we can transform the pion form factor Eq.Q into the compact parameter b space Issj]. 

F-niQ"^) = J[dXidh][dyjdh]^p{xi,h; iJ,f)f{xi,h;yj,h; ij.f)4j{yj,h; ^ x St{xi)St{yj) x 

exp{-S{xi,yj,Q,h,h;Hf)), (13) 

where fif = \ia{fif / Aqcd), [dxidh] = dxidx2d'^h5{l — xi — X2)/(167r^) and the hard kernel 



T{xi,h;yj,h; /i/) 



.e-^^-^^-^^-^^T{x,,k^i;yj,l^,;fif) 



(27r)2 (277)2 

The factor exp{—S{xi, yj, Q, b, h; fif)) contains the Sudakov logarithmic corrections and the 
renormalization group evolution effects of both the wave functions and the hard scattering 
amplitude, 



S{xi,yi,Q,h,h;fif) 



2 2 \ 1 U 1 

J2sixi,b,Q) + J2-^iyj^h,Q) - In^-— In 

Pi b Pi h 



/i/ 



(14) 



where b = ln(l/6AQC'j:')5 h = In^l/hAqcD) and s{x,b,Q) is the Sudakov exponent factor, 
whose explicit form up to next-to- leading log approximation can be found in Ref. Q|. In 
Eq.lfT^. St{xi) and St^yi) come from the threshold resummation effects and the exact form of 



each involves one parameter integration 
we take a simple parametrization proposed in Ref.jS, 

^ 2«T(3/2 + c 



35l |. In order to simplify the numerical calculations. 



x{l-x)Y 



(15) 



v^r(i + c) 

where the parameter c is determined around 0.3 for the pion case. 

To obtain the momentum projector for the pion, one may take the Fourier transformation 
of the bilocal operator matrix element defined in Eq.ljH]) jot, 

^/^f.z. N i_ N f^^^,^^ '^'aix^'^± ) „«^a(a;,kx) d \') 

6 



4 



J^75^7r(a;,k_L)-mg75 ^pp{x,'k±) - ia^ 



6 



9k I 



(16) 



where ip'^i^x, k±) = dipaix, k±)/dx. n = (1, 0, 0_l) and n = (0, 1, 0_l) are two unit vectors that 
point to the plus and the minus directions, respectively. Note we have used the parameter 
rriQ to replace the factor /i^ in Eq. (fTB|) . 

With the help of the above equations, the final formula for the pion form factor in the 
modified pQCD approach can be written as. 



16 



FAQ ) = -^^fnQ / dxdy bdbhdhasifif) X -ilj^{x,b; fif)ipl{y, h; ij,f) + 

y Jo Jo iz 



y 



P\2 



yijpix, b; ^ifWp{y, h; fif) + (1 + y)ipp{x, b; fif) 



3ijp{x,b;fif) 



6 



f{x,h;y,h]fj,f) x St{xi)St{yj) x 



exp{-S{xi, yj, Q, b, h; /i/)) 



(17) 



where x = {1 — x), y = {1 — y) and i/j'* {y , h; f) = d'ip*{y,h; fif)/dy. The first term in 
the square bracket gives the general twist-2 contribution and the remaining terms that 
are proportional to an overall factor ((^o)^/<5^) give the twist-3 contribution. The hard 
scattering amplitude T{x, b; y, h; fif) is given by 



(18) 



f (x, b; y, h-fif) = Ko (V^Qb) {e{b - h)K^ (^Qb) Jq {VVQh) + 
e{h-b)Ko{VVQh)lo{VVQb)), 



where the higher power suppressed terms such as (k^^/Q^) has been neglected in the nu- 
merator, Jo and Ko are the modified Bessel functions of the first kind and the second kind 
respectively. If taking out the threshold factors and absorbing the Sudakov factor into the 
definition of the wave functions, Eq.()17p agrees with Eq.(8) in Ref. 23[ (the factor before 
(^ijjp{x,b; fif)^jj*{y, h; fif)/Qj should be 3 other than 2 obtained there.). To ensure that the 
pQCD approach is really applicable, one has to specify carefully the renormalization scale 

lere are many equivalent ways to do so, a popular way 



fif in the strong coupling constant. T 
is to freez e a J Q'^) at lower Q^js, 3^, 
in Refs. , i-e. its value is taken as the largest renormalization scale associated with the 



37 



38, 



Here we take the scheme that is proposed 



exchanged virtual gluon in the longitudinal and transverse degrees. 



fj,f = ma.x{^/xyQ , 1/6, 1/h), 



(19) 



The Landau pole in the coupling constant at fif = Aqcd can be safely avoided in this way. 
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TABLE I: The full form of the LC wave function {x, k^) = ^^{x, k_L)x,r with the helicity function 
Xtt being included. ip-l'{x,'k±) stands for ^^1(2;, k_|_), ipf{x,'k±) and 'tpl{x,'k±), respectively. 



AiA 



1^2 



TT 



Ti 



V'(x,k_L) 



it 



ii 



^(x,kj 



ip{x,k_i 



Only the usual helicity components (Ai + A2 = 0) in the pion wave function have been 
considered in Eq. p7|l . From Eq. ()A2j) in the appendix, one may observe that the full form 
of the pion LC wave function have four helicity components (Table. HJ: namely, 



(20) 



By including the higher helicity components into the pion form factor, Eq. p7|) can be 
improved as 



16 

KiQ'^) = — tt/^Q^ / dxdy / bdbhdhasifif) x 
9 JO Jo 



|5:p(v^,^ai,a2)+ 
A1A2 



\ A1A2 A1A2 



2 E W;^,Ai,A2) 



f (x, b;y, h;/i/) x Stixi)Stiyj) x 



A1A2 

exp{-S{xi, yj, Q, b, h; /i/)) 



where t/'^' = d'lpl/dx and 



A1A2 



A1A2 



(21) 



A1A2 

A1A2 

In the above equation, because both the photon and the gluon are vector particles, the 
quark helicity is conserved at each vertex in the limit of vanishing quark masS|4^|. Hence 
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there is no hard-scattering amphtude with the quark's and the antiquark's hehcities being 
changed. For the hard scattering amphtude T{x, b; y, h; fij), we have imphcitly adopted the 
approximate relation for all the twist structures in Eq.(j2H), i.e. 

f{x, b; y, h; /i;)^^^^^ ^ _f^^^ b; y, h; (22) 



By ignoring the transverse momentum dependence in the quark propagator and applying 
the symmetries of the wave functions, especially the fact that k^) = ip^{x,li±), 

Ref. |82l| pointed out that the approximate relation Eq. ()22|l can be strictly satisfied. In 
fact, when the transverse momentum dependence in the quark propagator has been ig- 
nored, the Th depends only on one compact b-space, and Eq.(j22I) can be changed to a 
strict one, i.e. T{x,y,h; fify^~^^^ = —T{x,y,h;fify^^^K As is shown in Ref.|2], the trans- 
verse momentum dependence in the quark propagator will give about 15% correction at 
Q = 2GeV[2], so this effect can not be safely neglected. The hard scattering ampli- 
tude for the twist-2 contribution has been strictly calculated in Ref.^^ within the LC 
pQCD approach. One may find that when all the fcy dependence are included, strictly 
T{x,h;y,h; fify^~^^^ ^ —T{x,h;y,h.;fify^~^^^ and Eq.(j221) can be approximately satisfied. 
In the following discussions, we will keep the transverse momentum dependence in the hard 
scattering amplitude fully and use the approximate relation Eg . (|22j) to estimate all the 
helicity components' contributions to the pion form factor. 

Before doing numerical calculations, we would like to mention a few words on the value 
of rriQ. Based on the equation of motion of on-shell quarks, the authors used Ujt = m^/(m„ + 
mw) ~ 2.0GeV instead of rrtn for the twist-3 wave functions in Refs. Isl. l23l. 41 1 . A running 
behavior has been introduced in R,efs.|19l. l2flL 12 IL I22| and with this choice, one may find 



that the average value for /x^ over the intermediate energy region is around 2.5GeV. In 
Refs.Q Q a smaller phenomenological value ~ lAGeV, which is consistent with the 



result obtained from the chiral perturbation theory 



25, kj, is used to fit the B meson to 



the light meson transition form factors. Based on the moment calculation by applying the 



QCD sum rules, Ref. 2J] obtained rriQ = 1.30 ± O.OGGeV, which is very close to the above 
phenomenological value. So to be consistent with our model wave function constructed in 
the last section, we will take mg = l.SOGeV for our latter discussions. 

We show the twist-3 contribution to pion form factor Q'^F{Q'^) with all helicity compo- 
nents (i.e. using the full form of the LC wave functions ippi^y and ipl{x, k^)) calculated 

11 



0.35 



0.3 - 



0.25 




5 10 15 20 25 30 

Q^(GeV^) 

FIG. 2: Twist-3 contribution to the pion form factor Q-^F(Q^), where the second moment of 
ipp{x,'k±) or ^^(x,k_|_) is taken to be (^^) = 0.350. The dash-dot Hne and the dashed hne are the 
twist-3 contributions for the full form of the LC wave function with or without considering the fcy 
dependence in the quark propagator. The solid line is for the twist-3 contribution from the LC 
wave function that contains only the usual helicity component but is normalized to unity. 

within the modified pQCD approach in Fig.Q, where the second moment of k^) is 
taken to be (^^) = 0.350. One may observe that the transverse momentum dependence in 
the quark propagator will give about 25% correction at = 2GeV for the twist-3 contribu- 
tion, which is bigger than the case of the leading twist contribution. So it is more essential 
to keep the transverse momentum dependence fully into the hard scattering kernel for the 
twist-3 contribution. As a comparison, we also show the contribution from the twist-3 wave 
functions (i.e. k^) and k^)) that contain only the usual helicity component but 

are normalized to unity in Fig.Q- One may find the contribution from the twist-3 wave 
function that contains only the usual helicity component but is normalized to unity (the 
solid line) is larger than the contribution from the wave function with all the helicity com- 
ponents being considered (the dash-dot line). It is reasonable and is also the case of the 
twist-2 contribuion^^, because if one normalizes the valence Fock state to unity without 
including the higher helicity components, then the contribution from the valence state can 
be enhanced and become important inadequately. 
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IV. COMPARISON WITH OTHER MODELS FOR TWIST-3 WAVE FUNCTION 



As has been pointed out in Sec. Ill, the contribution from the twist-3 wave functions 
^/'p(x, k_|_) and k_|_) that contain only the usual helicity components but is normalized 
to unity is larger than the contribution from the wave functions k_|_) and k_|_) 

with all the helicity components being considered. However, as is shown in Fig.Q, since 
both of the twist-3 contributions have a similar behavior and are close to each other, the 
qualitative conclusions will be the same. And for easy comparing with the results in the 
literature, we will take the LC wave functions ipp{x, k_|_) and ipai^, k_|_) that only contain the 
usual helicity components for the discussions in the present section. 

Because of the end-point singularity, the twist-3 contribution depends heavily on the 
twist-3 wave function, especially on tpp{x,]<.±). In this section, we will do a comparative 
study on the twist-3 contribution from different type of ippi^y For this purpose, we take 
Eq.()17p to calculate the pion form factor, in which only the usual helicity components in 
the wave functions have been taken into consideration. 

The twist-2 and twist-3 wave functions tp-,^, ipp and ip^ may have different transverse 
momentum dependence, and for simplicity, we assume the same transverse momentum de- 
pendence for these space wave functions. For the transverse momentum dependence of the 
wave function, we take a simple Gaussian form, i.e. 

E x, kj_ = — - exp - , . , 23) 

g{x) \ 8f3^g{x) ) 

where A is the normalization factor, g[x) is either 1 or x(l — x). When g{x) = x{l — 
x), it is agree with the BHL prescription mentioned in Sec. II. After making the Fourier 
transformation, Eq. ()23j) can be transformed into the compact parameter b space as, 

^<^-''' ^ Mi'"">^ {-%^) M^^M . (24) 

where the upper limit (1/&) is necessary to insure the wave function to be "soft" 4J,|45|. 

Next, we consider the pion wave functions. The twist-2 wave function ip^i^, k_|_) with the 
prescription Eq.® can be written as 

V..(..k,)=^,exp(-^^gt^). (25) 
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TABLE II: The first three moments for the twist-2 and the twist-3 wave functions, where ah the 
full form of LC wave functions have the same BHL-like kx dependence, Tpl{iljl) = ■ip^(ip^)xTT and 





without Wigner rotation 


with Wigner rotation 










4'' 


4'' 








1 


1 


1 


1 


1 


1 


1 


ie) 


0.167 


0.350 


0.333 


0.391 


0.352 


0.176 


0.350 




0.060 


0.185 


0.200 


0.251 


0.197 


0.066 


0.185 



where the parameters can be determined by the normahzation condition of the wave 
function! 



^„.^|k,,^(,,k,) = l. (26) 

and some necessary constraints[31]. Taking the parameter values in Eq.ijH)), we obtain = 
1.187 X 10~^MeV~'^. The asymptotic form of twist-3 DA is the same as that of (pni^), 

and the end-point singularity coming from the hard scattering amplitude can also be cured. 
So for we also take its asymptotic form. For the twist-3 contribution to the pion form 



factor, the main difference for the existed results 3, 1201 121|, |23[ comes mainly from the 
different models for ipp{x, k^). The difference caused by the different model of ipaix, k^) (if 
all of them are asymptotic like) are quite small, so in the following, we will only consider the 
difference caused by different type of ipp{x, k_|_) and the contribution from k_|_) will be 
included as a default with the fixed asymptotic form for its DA and the same kx dependence 
as ipp{x, k±). 

In the asymptotic limit, (pp^ix) = 1, the end-point singularity coming from the hard 
scattering amplitude can not be cured, and the model dependence of (pp{x) is much more 
involved. Ta king the asymptotic DA and ignoring the fc^ dependence in the wave function, 



Refs 



.0 



21 



221 obtained a much larger contribution in a wide energy region 2GeV^ < 
< AOGeV"^, comparing with the twist-2 contribution. Using the model wave function of 
i^pi^x, k±) constructed in Sec. II, one may find that the twist-3 contributions are suppressed 
certainly. 

To study this effects more clearly, we compare our model with three different types of 
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k_L). In the literature, most of the calculations on the twist-3 contribution of the 
pion take (pp{x) as iljp{x,k±), i.e. without considering the intrinsic fcy dependence in the 
wave function, some examples for the electromagnetic pion form factor can be found in 



Refs. 



21, i3 



and examples for the B ^ tt 



However, as has been argued in several papers 
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brm factor can be found in Refs. 



23 



35 



42| 



the intrinsic transverse momentum 



dependence in the wave function is very important for the pion form factor and the results 
will be overestimated without including this effect. So in our comparison, the three different 
type of wave functions are constructed by adding a common simple Gaussian form fEq. ()23|) 
with g{x) = 1) to three different type of distribution amplitudes used in the literature^, i.e. 
the one of asymptotic behavior, the one in Ref . and the one in Ref . 24 1 respectively. 



4'\x,k^) = ^;exp('-^V (27) 



7/>f(a;,k^) = (l + 0.43C2'/'(2a;-l) + 0.09Cl/'(2x-l))A;exp(-^) , (2J 



V'f(x,k^) = (l + 0.137C2'/'(2x-l)-0.721Cl/'(2x-l))A;exp(^-^^. (29) 

The parameters and can be determined from the similar wave function normalization 
condition as Eq.^, = 7.025 x IQ-^MeV'-^ and P' = 168MeV. For the wave functions 
V^p'' — 1)2,3), the harmonic parameter /?' is different from that of '?/',r(a^, k^), however 
it can be taken as an effective /average value of the harmonic parameter with m = and 
g{x) = 1. The moments of the corresponding DAs are listed in TablelTTl 

We show the contributions to the pion form factor from the different model for ippix, k_|_) 
in Figs.(jnKISt'); where the contribution from our model wave function ipp{x, k_|_) with varying 
second moment (^^) is shown by a shaded band and the twist-2 contribution from k_|_) 
is included in Fig.(jnt) for comparison. Our present result for tjj^\x,'k±) (in dashed line) 
is much lower than the result shown in Ref. [2J|, since the value of /i^ used there has been 
changed to the present value of mg. One may observe that the twist-3 contribution is 
improved with our model wave function, and for the case of (^^) = 0.350, at about = 
SOGeV^, it is only about 45% comparing with the twist-2 contribution. This behavior is 



^ By using the prescription Eq.JHJ for the intrinsic kr dependence fEa. (|23|l with g{x) — x{l — x)), we can 
construct another three different model wave functions for ■ipp{x,\i±). However one may find that the 
moments of these three wave functions are too small and are out of the reasonable range obtained from 
the QCD sum rule, so we will not take them for our study. 
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Q'(GeV') 



Q'(GeV') 



(a) (b) 
FIG. 3: Different twist-3 wave function's contribution to the pion form factor, where the left 
diagram is for Q'^F{Q'^) and the right is for Q'^F{Q^). The dashed hne, the dash-dot hne and 
the dotted hne are for -i/'p^^ (x, k_|_), ^^^-'(a;, k_|_) and ^p^^(x,k_|_) respectively. The contribution from 
our model wave function '(/'^(x, k_|_) with varying second moment (.^■^) is shown by a shaded band, 
whose lower and upper edges correspond to (.^^) = 0.320, 0.370 respectively. For comparison, the 
twist-2 contribution from ■(/'^^(x, k^) is shown in solid line. 



quite different from the previous observations 
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20 



21 



22 



23|, where they concluded that 



the twist-3 contribution to the pion form factor is comparable or even larger than that of 
the leading twist in a wide intermediate energy region. 

As is shown in Figs.dnti'ISb), the twist-3 contribution from il)^\x^ k_|_) is comparable to our 
model wave function, which also has the right power behavior. We take a simple Gaussian 
behavior fEq. ()23|) with g{x) = 1) for the transverse momentum dependence in ip^^\x,'k±), 
i.e. a complete factorization between longitudinal and transverse momentum-dependence in 
the wave function. This Gaussian distribution behavior shows a strong dumping at large 
transverse distances, ~ exp(— 2/?^6^), while our model function with the prescription Eq.® 
has a slow-dumping with oscillatory behavior, ~ cos - x)b(3^ - 7r/4) /Vb. If we also 

take the simple transverse momentum behavior in our model wave function, i.e. the one as 
ip^^x, k_L), we find that the twist-3 contribution will be even lower, which is shown clearly 
by the dotted line in Figs.ijSti'jSb)- However as is shown in Fig.ijHb), "we can not achieve a 
right power behavior with k_|_), i.e. it drops down too quickly. 

Finally, with our model wave function for ipp{x, k), we discuss the model dependence of 
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the twist-3 contribution on the DA moments (^^"'). Here we take the second moment (^^), 
which gives the main contribution to (f)p{C,), as an example. Varying the second moment 
(^^) within a broader range, i.e. (^^) G (0.320,0.370), and adjusting the fourth moment 



(^^) to make <f)p{C,) has a closed behavior as the one that is obtained in Ref. 2J| (i.e. the 
dashed line in Fig.(P)), we can determine the corresponding parameters Ap, Bp and Cp in 
the wave function iljp{x,k.). The twist-3 contribution to the pion form factor with varying 
second moment {C,"^) has been shown by a shaded band in Figs.(jnKISb)- One may observe 
that the pionic twist-3 contribution increases with the increment of (^^) and all has a quite 
similar behavior on the variation of the energy scale Q^, i.e. as is shown in Fig.(jnb), the 
right asymptotic power behavior of order l/Q*^ has already been achieved at the present 
experimentally accessible energy region. 

V. SUMMARY AND DISCUSSION 

In this paper, we_have constructed a model wave function for ?/'p(x, k^) based on the 
moment calculation by using the QCD sum rule approach. It has a better end-point 
behavior than that of the asymptotic one and its moments are consistent with the QCD 
sum rule results. Although its moments are slightly different from that of the asymptotic 
DA, its better end-point behavior will cure the end-point singularity of the hard scattering 
amplitude and its contribution will not be overestimated at all. 

With this model wave function, by keeping the fc^ dependence in the wave function 
and taking the Sudakov effects and the threshold effects into account, we have carefully 
studied the twist-3 contributions to the pion form factor. Comparing the different models 
for ippix, k_|_), a detailed study on the twist-3 contribution to the pion form factor has been 
given within the modified pQCD approach. It has been shown that our model wave function 
ippi^, k_L) can give the right power behavior for the twist-3 contribution. With the present 
model wave function defined in Eq.® for iljp{x,'k±), our results predict that, at about 
~ lOGeV^, the twist-3 contribution begins to be less than the twist-2 contribution, and 
for the wave function '?/'p(a:;, k_|_) with (^^) = 0.350 at about = 30GeV^, it is only about 
45% comparing with the twist-2 contribution. This behavior is quite different from the 

where they concluded that the twist-3 contribution 



previous observations 



20 



21 



to the pion form factor is comparable or even larger than that of the leading twist in a wide 
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^ 0.4 


CD 



CM ^ 

o 



15 20 25 30 

Q^(GeV^) 



FIG. 4: Perturbative prediction for the pion form factor. The diamond hne, the dash-dot hne, 
the dashed hne and the sohd hne are for LO twist-2 contribution, the approximate NLO twist-2 



l3 



contribution 1 1(1. Il7l |. the twist-3 contribution _and the combined total hard contribution, respec- 
tively. The experimental data are taken from [4^. 

intermediate energy region up to AOGeV^. The higher helicity components (Ai + A2 = ±1) 
in the twist-3 wave function that come from the spin-space Wigner rotation have also been 
considered. The higher helicity components in the twist-3 wave function will do a further 
suppression to the contribution from the usual helicity components (Ai + A2 = 0), and at 
about = 5GeV^, it will give ~ 10% suppression. 

In Fig.Q, we show the combined hard contributions for the twist-2 and twist-3 contri- 
butions to the pion form factor, where the higher helicity components have been included 
in both the twist-2 and the twit-3 wave functions, and the twist-3 contribution has been 
calculated with our model wave function k^) with (^^) = 0.350. As has been pointed 
out in Refs. , the applicability of pQCD to the pion form factor can only be achieved 
at a momentum transfer bigger than ~ AGeV'^, so in Fig.(j3)), all the curves are started 
at = 4GeV'^. Together with the NLO corrections to the twist-2 contributions, which for 
the asymptotic DA, with the renormalization scale ur and the factorization scale /ij taken 
to be fi% = ^1^ = Q\ can roughly be taken as 10,|l7|, Q^F^^^ ^ [0.903GeV^)a^^{Q^), one 
may find that the combined total hard contribution do not exceed and will reach the present 
experimental data. There is still a room for the other power corrections, such as the higher 
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Fock states' contributions Q|; contributions etc.. Finally, we will conclude that 
there is no any problem with applying the pQCD theory including all power corrections to 
the exclusive processes at > a few GeV"^. 
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APPENDIX A: FULL FORM FOR THE LC WAVE FUNCTION 

By doing the spin-space Wigner rotation, we can transform the ordinary equal-time 
(instant-form) spin-space wave function in the rest frame into that in the LC dynamics. 
After doing the Wigner rotation, the covariant form for the pion helicity functions can be 
written as 0,0], 

Xn{x, kx) = U{pi, Xl)j5v{p2, A2), (Al) 

where pi = (x, k_|_) and p2 = (x, — k_|_) {x = 1 — x) are the momenta of the two constituent 
quarks in the pion, Mq = (m^ + k'^)/x{l — x) and the LC spinors u and v have the Wigner 
rotation built into them. Then the full form of the LC wave function can be written as 

t/j^ix, kx) = tPix, k^)xAx, kx), (A2) 

where the momentum space wave function ip{x,k±) represents ■iljp{x,'k±), 'iljT,{x,k±) and 
tpcrixjkx) respectively. Because all the LC wave functions can be dealt with in a similar 
way, here we only take ipp{x,]<.±) that is defined in Eq.(jHI) as an explicit example to show 
how to determine the parameters in the full form. 

The full form of LC wave function ipjjix, k^) contains all the helicity components' contri- 
butions and its four components can be found in Table|l] The parameter values built in the 
wave function ip^^x, k±) can be done in a similar way as for the wave function of ippi^y k±) 
that contains only the usual helicity components, i.e. 

v4p = 4.088 X 10~^Me■l/-^ 5^ = 1.077, = -4.317 x 10"^ (A3) 
m = 309.6MeV, /? = 395.9MeV, for (k^^) ^ {367MeVy (A4) 
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where the parameters m, (3 are determined by the wave function normahzation condition 



and some necessary constraints 



3l|, and the values of Ap, Bp and Cp are determined by 



requiring the first three moments of its DA to be the values shown in Eq. ()ll|) . From the 
wave function Eq. ()A2|) . we obtain 



ApTTlP 



TX--' 



3/^2(1 -a 



;i + BpCl'\i) + CpCl'\i)) exp 1 - Erf 



(A5) 

where the error function Erf{x) is defined as Erf{x) = Jq e~^^dt. One may find that 
4>p{C,) is almost coincide with (f)p{C,) that is shown in Eq. ()l()|l . and for simplicity, we can take 
the approximate relation, 4>^{C,) ~ 4>p{0- It is reasonable because we have adjusted the 
parameters in both DAs to have the same moments and due the fact that the momentum 
space wave function ipp{x, k^) is an even function of k_L, one may find that the higher helicity 
components in ipjii^: do not contribute to (p^iC)- 
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